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I. INTRODUCTION 

Several years ago, Vafa and Witten[l] (hereafter referred to as VW) provided a simple, elegant argument for the 
absence of spontaneous parity violation in theories containing only vector-like fermions; this result, however, has been 
criticized on various grounds dealing with mathematical niceties underlying the proof. It was argued in [2] and later 
in [3] that the expression for the free energy on which the VW argument is based has an ill-defined infinite- volume limit 
and is therefore insufficient to determine the absence of spontaneous parity violation. Previously, Aoki[4] had noted 
the possibility that lattice QCD has a parity- violating phase which disappears in the continuum limit; this apparently 
presents a problem for the VW result since the original proof can be applied to the lattice-regulated theory as well as 
to the continuum. Hence the Aoki phase, even if only present at non-zero lattice spacing, would directly contradict 
the results of [1] provided the assumptions made by VW are satisfied by the lattice regulated theory at finite spacing. 
This model was studied carefully by Sharpe and Singleton[5] who argued that the Aoki phase is not ruled out by the 
results of VW since their arguments appear not to be applicable for the case of fermionic order parameters. Finally, 
reference [6] raised the possibility of parity- violating operators that evade the proof of VW for systems in a thermal 
bath. In this paper, we will re-examine these complaints and find that, for physically relevant systems, they do not 
invalidate the results of VW.[l]In a related paper, [7] we argue that the proof, even under the assumptions stated, 
does not however exclude spontaneous parity violation under certain circumstances. 

We organize this paper as follows. In section II we review the original VW argument. In section III we examine the 
large volume limit of the theory and show that, as claimed in [2], it is well defined only for the case where spontaneous 
symmetry breaking docs not occur; yet we will argue that this deficiency does not invalidate the results of VW. 
We revisit the lattice arguments in section IV where we argue that the Aoki phases (at least for a large number of 
fermions) do not represent absolute minima of the effective potential but only metastable states 1 . In section V we 
consider the use of fermionic operators as order parameters for a possible parity- violating phase transition; in contrast 
to the results of [5] we find that this type of order parameter exhibits the same behavior claimed by VW. We then 
argue that there is a simple and useful fermionic operator that can be better used as an order parameter since it 
preserves chirality (the operators used in [4, 5] violate both parity and chiral symmetry). Parity violation at finite 
temperature is briefly revisited in Sect. VI, and parting comments are presented in section VII. A mathematical result 
is relegated to the Appendix. 



This should not be used to completely dismiss these states as irrelevant: recent experiments in heavy-ion collisions create finite-volume 
systems that may relax temporarily into parity- violating configurations [8] 
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II. THE VAFA-WITTEN ARGUMENT 



In this section we summarize the original VW argument. The idea is to consider a theory with action S that is even 
under parity transformations and contains only vector-like fcrmions (so that the corresponding fcrmionic determinant 
is real), and to add a parity- violating term X so that the full Minkowski-space action becomes 



Stot — S + XX, 



(1) 



where A is a small real parameter and X is Hermitian. We will assume (without loss of generality) that X is parity 
odd; all gauge-fixing terms are assumed to be contained in S. The bosonic terms in S are assumed to be real (this 
is always possible for all linear gauges) 2 , (see for example [9]). We will concentrate on theories that do not contain 
scalar fields. 

We first consider the case where X is a purely bosonic operator. Some comments concerning scalar theories are 
provided in section VII. We first obtain the free energy by Wick-rotating S in (1). While not without ambiguities, to 
be precise, we will follow the procedure of [10]. Any parity- violating operator contains an e tensor and this implies 
that the corresponding term acquires an additional factor of i upon Wick rotation, thus we find 



iS tl 



St 



iXX, 



(2) 



with Se and Xe real (the subscript E will be used to denote the quantities in Euclidean space). The "free energy" 
W is then 



-w(\) 



'dA]e 



Det(0 E + AI), 



(3) 



where A denotes the gauge fields, and S E oson the purely bosonic contributions to Se- 3 We take the mass matrix 
M diagonal with non- negative eigenvalues (as we are not interested in explicit P violation); the Euclidean-space 
fermionic covariant derivative is denoted by De which is anti-Hermitian in our conventions [10] (the gamma matrices 
are Hermitian). The fermionic determinant is then a real functional of the gauge fields for vector-like couplings [11]. 
Consider now the case where X is a fermionic bilinear 4 , 



X = tpOip. 

After a Wick rotation and to first order in A the fcrmionic integration yields 



-W(X) 



;dA]e- s B° aon+lAX s lf) Det(^ B + M); X 



(off) 



-iTr 



Pi 



M 



:0, 



(4) 



(5) 



where Oe = 7°7 5 07 5 7° 5 . A naive calculation shows that X E is real if O e = —^Oe^ 5 - In particular note that 
this holds for the simple cases Oe = «7s and Oe = PeI^ ■ There is, of course, a possibility that these rather formal 
results are flawed by the mathematical subtleties involved in properly defining xj^ ff '[5]; we discuss this possibility it 
in section V. 

With these preliminaries the proof is based on the claim [1] that all parity-violating operators correspond to real 
Xe or Xe , in which case 



W(A) > W(0), 



(0) 



as a result of Schwartz's inequality. The validity of this result then rests on the requirement that the measure is 
positive definite, that is, Det(pE + M) > 0. This is true for the case of vector-like fermions [13], but can be extended 
to some more general types of fermionic Lagrangians; we return to this point in section IV. Note that in proving (6) 



2 For this it is convenient not to introduce ghosts and deal with the gauge-fixing determinant directly. 

3 For constant A, W-^A) will be proportional to the space-time volume. 

4 When X contains products of bilinears these can be turned into a sum of bilinears through the introduction of appropriate auxiliary 
fields by using the Hubbard-Stratonovich trick [12]; the arguments below carry through provided these auxiliary fields are vectors; this 
will be investigated in detail in [7] . 

5 It is understood that Oe is obtained by replacing Do — > —iD^, 7 — > 7^, ■y k — > — 17^ in O. In our conventions yjjj< = +7^ and 
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we used the fact that W is real since X is odd under parity and the integration involves both a configuration and its 
parity conjugate. 

Since any explicit violation of the symmetry (in this case parity) results in an increased free energy, the lowest 
energy ground state will be parity symmetric. To see that eq. (6) implies vanishing of the expectation value of X , X, 
it is useful to pass to the effective potential V(X) 

— — — dW 

V(X) = W(X)-XX, withX=— . (7) 

OA 

Inversely, using A = —V'(X), we can write this as 

W(X) = V(X) + XX = V(X) - XV (X). (8) 

We want to show that, in the limit A^O*, X^O. Suppose, on the contrary, that 6 as A — > + , X = v ^ 0. This is 
necessarily at a minimum of V(X). 7 If we adopt the usual convention that W(0) = 0, then V(v) = as well. The 
inequality eq. (6) therefore implies that 

V(X)-XV'(X) = {V(X) - (X-v)V'(X)} - vV'(X) > 0. (9) 

Since v is a minimum, V (v) = 0. If we expand in a neighborhood of X = v, the term in the curly brackets contains 
at least one power of (X — v) greater than the last term and, moreover, since v is a minimum, the leading term in 
the expansion of V'(X) will be an odd power of (X — v). (For V"(v) > 0, it will be linear in (X — v).) But then the 
inequality will be violated in a neighborhood of X = v. Therefore, the only possibility is that v = 0. 8 . 

The same conclusion can be obtained starting form the Hamiltonian derived form (1) 9 . Indeed, if the ground state 
is not parity symmetric, then, when A = 0, there will be two degenerate lowest-energy states |±) that are parity 
conjugate of each other. Now, the Hamiltonian corresponding to (1) is of the form H + AS where H is parity even 
and H parity odd; a straightforward application of degenerate perturbation theory shows that the term AS lifts the 
degeneracy, so that the lowest energy state for the system equals 

E a -\Xe\+0{X 2 ); E = (±\H |±>; e = (+|S|-) (10) 

that decreases as a function of A. This contradicts (6) and implies that there are no parity-degenerate vacua: the 
symmetry is not broken. 



III. LARGE VOLUME EFFECTS 



The results summarized in the previous section were criticized in [2] based on the following observation. Suppose 
we construct a Landau-Ginzburg functional of the order parameter defined as Xg/V, where V denotes the spatial 
volume of the system. Explicitly, 



-W 



dp J [dA]e- s ^° n+lXXE Det(# B + M)S(p - X E /V), 

dpe -HP)+iVX P) 



(11) 



where T is presumably proportional to V and we assume it can be expanded in a power series in p. Since the original 
action is parity symmetric T will be even in p. For definiteness we will consider the case 



Hp) = v 



(12) 



6 The case A — » is similar. 

7 Since the action Se is parity symmetric, and X is parity odd, W(A) = — W(— A), and V(X) = V{— X). So, —v will also be a minimum, 
but we simply pick one of them. 

8 The same result is obtained when one defines the effective potential in the region |X| < \v\ using the Maxwell construction [15]. 

9 The following is a clarification of an argument presented in [1] 
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hence 

e~ w = (^6)- 1 / 4 X(g) 



/OO 
dx exp 
-oo 



1 4 2/^ 2 



(13) 



where 

g = A(y 3 A>) 1/4 ; a=^a(A 2 6)- i/J . (14) 
The integral X is evaluated in the Appendix, we find, for small |A|, 

X ~ ze -^[a 2 /(4fc)+A 2 /(4a) + ...] cos | A | ^R7^ + </>) 5 |A| « 1, O < 

~ ze -via 2 /W + \MV^~ b +^ 2 /^)+-l j |A| « 1, a > (15) 

where z, </> are constants. 

For the case of no spontaneous symmetry breaking (a > 0) the free energy F is well defined and satisfies (6) though 
it is not analytic in A. When the symmetry is spontaneously broken the free energy suffers from the ailments described 
in [2]. Note that the bound (6) is satisfied in either case. 

Since the free energy is ill defined when a < (though (6) is still obeyed), the VW argument is incomplete in this 
case. The problem is that the partition function oscillates, corresponding to a complex W; we now argue that this 
will not occur in any physical system. Indeed, a reasonable requirement for a system is that the Hamiltonian be an 
additively renormalizable operator that is bounded from below. In this case exp(— W) equals the (3 — » limit of the 
partition function Z = Tr e~@ H . For for any finite /3, the standard expression for Z in terms of energy eigenvalues 
shows that it cannot be negative, provided we assume that the exponential decrease in the contributions from large 
(renormalized) energy states is sufficient to guarantee convergence of the trace. Systems obeying these restrictions 
will have Z > 0, which then implies exp(— W) > 0. In such cases the possibility of negative a (or any set of parameters 
leading to a complex W), where parity is spontaneously broken, will not be realized. 

This same issue was discussed in [3], where it is claimed that W(A), when defined appropriately, has the same 
dependence on A independent of the sign of a. Although we find that this is not the case, the claims of this reference 
would be validated by (15) provided the free energy is a concave function of A, as is indeed the case when a > 0. 

This, however, is not the case a < 0. If we adopt the prescription of [3] of defining the free energy through an 
average over the volume, the results are ambiguous. For example the free energy per unit volume, 

/ = - Jim ilnl (16) 

V — >oo V 

would equal / = a 2 /(4b) + X 2 /(4a) + • • • provided we define 

lim = (17) 

for any fixed real constant c. In this case the conclusions of [3] follow: we cannot use the A dependence of the 
free energy to determine the presence or absence of symmetry breaking. However this is an ad-hoc choice: consider 
V„ = (n + 1/2)tt — e~ U7rn for a real constant u and a positive integer n. Then letting V = V n and taking the infinite 
volume limit by letting n — > oo, we find 

r IncosK . 

lim — i — u, (18) 

having taken the branch of the logarithm such that lncosy„ = mr(i — u). With these choices the limit is certainly 
non-zero. The point is that the large volume limit of (lnX)/V is ill defined. 10 It is possible to define the limit so that 
the troublesome oscillatory term does not contribute, but it is unclear whether this is the physically correct definition 
of the free energy. 

We conclude that the concerns raised in [2, 3] are valid in that the partition function would not be positive definite in 
a system exhibiting spontaneous parity violation. However, this situation will not arise for systems with a Hamiltonian 
bounded from below such that the density of states does not grow faster than the exponential of the energy. 



For a different choice of the branch cuts the imaginary term would be absent. 
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IV. LATTICE CONSIDERATIONS 



In a series of publications [4] Aoki proposed a phase diagram for lattice QCD that is markedly different from the 
standard case (both quantitatively and qualitatively) for non-zero lattice spacing (all such peculiar effects disappear 
in the continuum limit). In this scenario the theory has two phases, one where the 3 pions are massless, being the 
Goldstone bosons produced by the spontaneous breaking of the usual chiral symmetry. In the second phase flavor 
symmetry is believed to be spontaneously broken down to U(l) generating two Goldstone bosons (for the case of 2 
massless quarks) so that the third pion is massive in this phase. The transitions between phases are of the second 
order so that the third pion mass vanishes smoothly as the phase boundary is approached. In this second phase parity 
is also spontaneously broken leading to a non zero expectation value of the order parameter (tpij^i/j) . 

Though Aoki's original arguments were made within the lattice-regulated theory, the results of Sect. II would 
be equally applicable to this case[5] provided the various assumptions made are satisfied. Should this be the case, 
the presence of a parity- violating phase at finite lattice spacing casts doubt on the results of [1], even though the 
parity-asymmetric phase disappears in the continuum. 

The lattice arguments in ref. [4] are supported by two large-N calculations, one for the 2-dimensional Gross-Neveu 
model[I6] and one for 4-dimensional QCD[I7]; both are based on finding minima for the effective potential that do 
not respect parity. In addition there is also some numerical evidence based on the Monte-Carlo evaluation of several 
observables. We will revisit these calculations in the paragraphs below. 



A. The lattice Gross-Neveu model 



This model contains N spinor fields ipi where the doubling problem is fixed using Wilson fermions. The dynamics 
is defined by the lattice action (we take a unit lattice spacing) 

N N 

S = 2 J2iH(nj { [ftfc] (n) + \r [aPi] (n) + [a(n) - m{n)^} ^(n) } + ^ E H 71 ? + < n ?} > ( 19 ) 

i=l ^ ) ^9 „ 

where g is iV-independent and the derivatives denote the usual lattice operations[19]. In the large N limit and for 
constant a and n, the fermionic determinant yields a function which has a minimum for a non-zero value of n. Even 
though quantum fluctuations are known to destroy this non-trivial vacuum[20], this calculation does suggest the 
possibility of the richer phase diagram described in [4] . 

These results for the Gross-Neveu model do not contradict general argument presented in [1] . This is because this 
model does not satisfy an important assumption made in deriving the VW result: since the fermions do not couple 
vectorially, the fermionic determinant associated with (19) is not positive definite. While this model can be used as 
guidance in obtaining the phase structure of QCD on the lattice, it does not provide a violation of the VW results. 



B. Lattice QCD in 4 dimensions at large N 

Another example considered by Aoki is the case of the large N limit (N again denotes the number of fermions) 
in lattice QCD in 4 dimensions using Wilson's formalism. Denote by cr(n) = ip(n) <g> ip(n), then, following [4], we 
consider the effective potential for these operators in the case where they are position-independent in the large N and 
strong coupling limits. Replacing then 



the effective action for this order parameter is (we take again unit lattice spacing) [17] 

1 + yjl -4cr 2 sin(6>) 2 



S eS = ANfl { ma cos(0) - ln(cr) - 2 



y/l-Aa 2 sin(6>) 2 - 1 - In 



(20) 



(21) 



where denotes the lattice volume, m the bare fermion mass, and we have chosen the Wilson fermion parameter 
r = 1. The effective potential is minus the above expression. We will also replace Nina = (N/2) In a 2 and avoid a 
spurious imaginary contribution to S e s. This expression has a series of local extrema at the points 



(i) : cos0 = ±l, <j = ±l/m. (ii) : cos9 = ±m 



16 — m 2 ' 



± 



16 — m 2 



(22) 
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Points (ii) are local minima of S e g for m 2 < 4 but arc saddle points for m > 2; (ii) are saddle points for m 2 < 16 and 
unphysical otherwise. It is remarkable that the curvature of the effective potential at points (ii) is positive in the 9 
direction, corresponding to a positive pion mass parameter. 

None of these points, however, represent the global minimum of the effective action. In fact we have the behavior, 

Seff — ► +oo, a — > 0, 

— > — co, a — > oo, 9 = n. (23) 

Hence points (i) represent meta-stable states where parity is broken, the stable vacuum at a = is parity even. 
Note that the stability of the a = vacuum is determined solely by the term In a and may be altered by finite N 
corrections. 

This still does not explain why the numerical evaluation of (V'75'0) yields a non-zero result. This is understood by 
the observation [4] that the Wilson action does not lead to a positive fermionic determinant, which follows form the 
fact that the Wilson term, does not anticommute with 75. It is possible to generalize somewhat the results of [13], 
and to show that the Wilson action leads to a positive fermionic determinant, but only when the fermions carry a 
real representation of the gauge group [7], which is not the case for QCD. Hence the presence of spontaneous parity 
violation in QCD with Wilson fermions does not invalidate VW for finite lattice spacing, as the theory does not 
comply with the assumptions made in [1]. On the other hand, the lattice results do indicate that (ipfsip) —* in the 
continuum limit, as demanded by VW. In this sense the results in [4] support those of [1]. 



V. FERMIONIC ORDER PARAMETERS 



In a calculation related to Aoki's, Sharpc and Singleton[5] considered the case of QCD with two flavors to which 
they added a source for a parity-violating order-parameter. Then they studied the vacuum expectation value of this 
order parameter as the source strength was set to zero. The fermionic part of the Lagrangian in Euclidean space is 

£fcrm = *(fe+m + te- ! ^ T3 )* ! (24) 

where ^ denotes the fermion iso-doublet (assumed degenerate in mass) , and h the source strength. The Euclidean co- 
variant derivative operator^)^; is anti-Hcrmitian (in our conventions the Euclidean gamma matrices are Hcrmitian[10]), 
and the Pauli matrix T3 acts on the isospin variables. 

The possibility of spontaneous parity violation can be studied considering the expectation values of operators of 
the form 

( r a) =lFT^ a , a = 1,2; (25) 

where a labels the isospin components and T is a linear combination of the unit matrix, 175 , T3 and ^375 . To obtain 
the general form of the expectation value, 

= ^ { T Pe + m + hcxpl-iO-yMaa] } ' (26) 

(no sum over a) we follow [5, 21] and consider the above system in a finite volume; we then expand the above trace 
using the eigenfunctions of Pe which we denote by ip n . Explicitly, 

Pe^Pu = iK<Pn, (27) 

where we used the fact that the Euclidean Dirac operator is anti-Hermitian. As usual [22] we will treat the zero modes 
separately. 

When A„ 7^ the mode (p n is paired with 75<p rl for which 

Pe (75<^n) = i(-A„) (75¥>n) , (28) 

hence in the {tp n , 75 ^n} subspace we have 

P E + m + he^^ (iXn + m + hcose T^sintf \ 

r \ ^fihsmU — i\ n + m + h cos 9 J v ' 

and we can replace 75 — > t± (denoting the usual Pauli matrix). Then the corresponding contributions to the trace 
(26) are 

f 1 1 2(hcos9 + m) J 1 1 2/isin6> 

1 1 Pe + m + /ie=F«75 J A 2 + A 2 ; \ (l7s) p E + m + /ie=F«7* J ~* T A 2 + A 2 5 ( } 
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where 

A 2 = h 2 + m 2 + 2hm cos 9. (31) 
For the zero modes A„ = we can choose <p n to be an eigenstate of 75, 

75V« = Xn(fin, xl = 1, (32) 

In this one-dimensional subspace we can replace 75 — > Xn an d 

+ m + he Tl0 ' t5 -> /le^"*" + to. (33) 
Then the corresponding contributions to (26) become 

^l ^ + rn + ^^ J"" A 2 5 ^l^ ^ + m + fte^ /^'X" A 2 " (34) 



Adding all terms we find 



1 1 v-~* h cos 9 + to /i sin 

TH 1- : 77-^ = > —tt: T77-±W- 



£ + to + /ie=F^75 f A 2 + A 2 A 2 ' 

' n ' 

{1 ^-^ /i sin 9 + m /z cos # + m 

( * 75 W^ + fe^ J = A 2 +A 2 +^ A 2 ' (35) 

where the sum is over all modes (with positive negative and zero eigenvalues), and v denotes the index of Pe, 



E 



where this sum is over the zero modes only. 

Using these expressions we find, for the operators in (24), 



Xn, (36) 



,(a)\ h cos 9 + to v A .. . v hsu\9 
' = Aft A 2 + A 2 _ aa ?2 A 2 ' 

n n 

( a )\ /isinfl^ A .vhcos9 + m 



AO ^ A 2 + A 2 n A 2 
(no sum over a) where fi denotes the space-time volume. In particular, 



n . hcos9 + m 

drx ^75* ) = 2iv — (38) 



Taking the infinite volume limit and then letting h approach zero, we obtain 



(fx *i7 5 * ) = 2i— 

TO 



h=0 

i 

8^ 



d xTr 



(39) 



The authors of [5] obtain zero instead of the above expression since they did not include the contribution from the 
zero modes explicitly. Note also that both isospin components contribute equal amounts to the right-hand side of 
(39). The fact that we obtain a non- vanishing expression in (39) docs not imply that parity is broken: this will be 
true only if this non-zero value survives the integration over the gauge fields. That integral cannot be evaluated, even 
formally; yet the fact that the right hand side of (39) is purely imaginary allows us to use the argument leading to 
(6) which implies that, in fact, parity is not broken. 

It is worth noting that (39) is consistent with the anomaly equation in Euclidean space provided we assume that 
there are no massless excitations, so that the volume integral of the divergence of the axial current vanishes. We also 
note that the corresponding expression for the parity-even bilinear gives 

/9m 
^PAW^—; (40) 
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where pa denotes the spectral density per unit four- volume for the operator Pe- This expression reproduces the 
corresponding result of [5]. 

If we add to the Lagrangian a term 9Ti F F / (4tt) 2 (where F denote the field strength of the gauge field in the fermionic 
covariant derivative and F the corresponding dual) then (39) can be interpreted as stating that the expectation value 
on the left-hand side can be compensated by an appropriate shift in 9. For the case of low-energy QCD this is precisely 
the same result obtained using a chiral Lagrangian: a constant shift in the rf field can be compensated by a shift in 
9 [23]. 

The expression (39) validates the claims of [1] in that the expectation value is purely imaginary. We also note that 
(39) is ill defined in the massless case where it has a singularity. This problem is connected to the fact that the 
operators under consideration, o\^ r violate both parity and chiral symmetry. We can repeat the calculations for an 
operator that violates only parity such as 

=Wi lb pt>% (41) 

In this case we find 

(0 { £>) = -i(T 3 ) aa hsm9 J dXp A (X) x2 2 ^ m2 , (42) 

that vanishes as h — ► 0, but also vanishes when 9 = when the theory is parity symmetric. For this operator the 
expectation value is zero at the level of the fermion integral, in contrast with (39). 

It is worth pointing out that Op is not related to 0^ in (25) through the use of the equations of motion, in 

fact, if we replace p —>p + im in (41) Op remains unchanged. This shows that this operator can be eliminated by 
performing a chiral rotation on the fermion fields, and is then equivalent to a purely bosonic operator proportional 
to the index of the Dirac operator Pe [24]. 

The authors of [5] also consider the behavior of the lattice version of this theory close to the continuum. They find 
that the leading contributions to the effective potential for the meson field £ = a + ir ■ tt, (a 2 + tv 2 = 1), equals 

Kff(S) = -cia + c 2( 7 2 , (43) 

with c 2 approximately constant (independent of m) while C\ is expected to be a linear function of m. The crucial 
issue here is the sign of c 2 : if c 2 < the potential has a minimum at a 2 = 1 that corresponds to a parity-symmetric 
vacuum where the pion field gets no expectation value. If, on the other hand, c 2 > parity will be broken 11 . In 
view of the above discussion that supports the results of Ref. [1], and suggests that the parity- violating vacua of the 
large-N, strong coupling lattice theory are not absolutely stable, we believe that, in fact, c 2 < as well. 

VI. PARITY VIOLATION AT FINITE TEMPERATURE 

In [6] it is noted that systems in a heat bath have an available additional time-like vector (the temperature 
vector), that can be used to construct P- violating operators that apparently do not acquire the crucial factor of i 
upon a Wick rotation. An example is the operator e a ^ s ngTi \D a Fp v n v F llL n ll, \. 

In order to examine this claim we consider the case of time-dependent field theory at finite temperature. In this 
case it is well known[25] that in the functional integral exponent, the time integration should be carried along a 
complex path starting at — oo and ending at — oo — if3. The temperature vector is, in the rest frame for the system, 
the tangent to this path, and so it is, in general, complex. For the purposes of evaluating the expectation value of an 
order parameter the path can be taken parallel to the imaginary axis and, as a result, is purely imaginary. Since 
all operators in the class mentioned above are odd in n^, they do receive a factor of i and will not violate the claims 
of [1]. 

VII. COMMENTS AND CONCLUSIONS 

The above calculations provide a re-examination of the original Witten-Vafa result concerning the absence of 
spontaneous parity violation in theories with only vector-like fermion couplings. We also reviewed the various criticisms 
of that result and concluded that none of them is strong enough to invalidate the conclusions of [1]. 



11 In addition, the phase structure suggested by Aoki is reproduced by this model when ci > 0. 
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Theories with scalars can evade the VW result in two ways. Firstly some pseudoscalar fields (if present) can acquire 
an expectation value, but such order parameters do not receive a factor of i when Wick rotated. Secondly, in such 
theories the fermionic determinant is not necessarily positive definite (e.g. in theories containing solitons.[26]). 

In confining and parity-conserving theories containing fermions and gauge bosons, the low-energy excitations are 
often scalars and fermions (as is the case for QCD). The comments above suggest that the low-energy effective theories 
for such models can, in principle, violate parity spontaneously. In fact the standard chiral non-linear sigma model[23], 
with the addition of an electromagnetic mass [27], allows parity- violating vacua for a certain range of parameters. On 
the other hand, the VW argument implies that there is no parity violation in the underlying theory Of course, one 
could simply argue that the VW result for the underlying theory merely requires that the low-energy parameters 
are such that no parity violation occurs. Still, this situation opens the possibility that there might be some non- 
perturbative effects that allow such theories to evade the results of VW. An investigation of this possibility lies 
beyond the scope of the present publication, though we intend to discuss it fully in the near future [7]. 

Finally, we note that, as briefly mentioned in [1], the same arguments cannot be used to rule out spontaneous 
CP violation in theories without scalars as there are CP-violating operators that remain real upon Wick rotation 
(e.g. F^F^pFp^fabc where F denotes the non-Abelian field strength, and / the group structure constants). So, even 
though it is known that such theories cannot violate CP explicitly, [28] the possibility of their exhibiting spontaneous 
CP violation remains. 
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VIII. APPENDIX 



In this appendix we evaluate the integral 



J —i 



dx exp 



aq 



2 / 3 x 2 +i 



q:r 



(44) 



in the limit of large \a\. Our strategy will be to construct a differential equation for X whose solutions can be obtained 
by a method similar to the WKB approach. The disadvantage of this approach is that the differential equation has 
solutions that do not correspond to X. In order to extract the relevant ones we will use the following property, 



HI < 



exp 



1 



aq 2 ' 3 x 2 



yJ\a\qWK 1/4 (W /3 /2 



2 1 /V/3 



that follows from the Schwartz inequality. 

Integrating by parts, we find (a prime denotes a q derivative) 



uX" + vX' + wX = 



where 



-a 2 q 1 ' 3 + 



8o 



6aq - (27 - 32a 3 )< 7 7 /3 



-a (4a 3 - 9) q 2 ' 3 - -a 2 q- 2 ' 3 - 



(224/9)a 3 9 



3„l/3 



^V 5/3 - ^<r 1/3 + |k<r 1/3 



6aq - (27 - Z2a 3 )q 7 ' 3 

16a 3 ,\ 32 a 2 (aq- 2 ^ + iq 2 ' 3 ) 



27 



1 3 



9 6aq - (27 - 32a 3 )g 7 / 3 



In the large q limit the solutions for X can be obtained by using the Ansatz 

J = exp (-yq i/3 ) 



(45) 



(46) 



(47) 



(48) 
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that solves (46) for large q provided 

64y 3 - 128a 2 ?/ 2 + 16a(8a 3 - 9)y + (16a 3 - 27) = (49) 

Of the 3 solutions to this equation one is real for all a; we call it yo. The two others are complex conjugate of each 
other for a < 0.944941 and real for a > 0.94491; we denote them by y± with the sign associated with the sign of the 
imaginary part. The asymptotic behavior of these solutions is given in Table I. 



solution 


a — ► +oo 


a — > — oo 


yo 


a 2 + v / 2^+ 1/(8q) H 


-l/(4a) + — 


y+ 


-l/(4a) + --- 


a 2 +iV2R+ 1 /(8a) + --- 


y- 


q 2 + V2a + l/(8a) H 


a 2 -i N /2H + l/(8a) + --- 



TABLE I: Asymptotic behavior of the solutions to the differential equation (46) 

The required solutions should satisfy (45) which implies that the corresponding y should have a positive real part. 
In this way we find, for large |a| 

1 ~ Ze -7 4/3 [« 2 +i/(M+-] cos (y2H<7 4 / 3 + A; a -> -oo 

X ~ Ze - g 4/3 [a 2 +v^+i/(8 Q )+...]. a ^ +co ( 50 ) 
where Z, (/) are constants. This is the result used in the text. 
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